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Abstract. We present cosmological parameter constraints on flat cosmologies dominated by dark energy using various cosmo- 
logical data including the recent Archeops angular power spectrum measurements. A likelihood analysis of the existing Cosmic 
Microwave Background data shows that the presence of dark energy is not requested, in the absence of further prior This comes 
from the fact that there exist degeneracies among the various cosmological parameters constrained by the Cosmic Microwave 
Background. We found that there is a degeneracy in a combination of the Hubble parameter Hg and of the dark energy equation 
of state parameter Wq, but that Wq is not correlated with the primordial index n of scalar fluctuations and the baryon content 
flb'i^- Preferred primordial index is n = 0.95 ± 0.05(68%) and baryon content Q.bh- = 0.021 ± 0.003. Adding constraint on 
the amplitude of matter fluctuations on small scales, erg, obtained from clusters abundance or weak lensing data may allow to 
break the degenaracies, although present-day systematics uncertainties do not allow firm conclusions yet. The further addition 
of the Hubble Space Telescope measurements of the local distance scale and of the high redshift supemovae data allows to 
obtain tight constraints. When these constraints are combined together we find that the amount of dark energy is O.T^Jjg" (95% 
C.L.) and that its equation of state is very close to those of the vacuum: Wq < -0.75 (> 95% C.L.). In no case do we find that 
quintessence is prefered over the classical cosmological constant, although robust data on erg might rapidly bring light on this 
important issue. 
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1. Introduction 

The determination of cosmological parameters has always been 
a central question in cosmology. In this respect the measure- 
ments of the Cosmological Microwave Background (CMB) 
anisotropies on degree angular scales has brought one of the 
most spectacular results in the field: the flatness of the spa- 
tial geometry of the universe, implying that its density is close 
to the critical density. Although, during the last twenty years 
the evidence for the existence of non-baryonic dark matter has 
strongly gained in robustness, observations clearly favor a rel- 
atively low matter content somewhere between 20 and 50% of 
the critical density, indicating that the dominant form of the 
density of the universe is an unclustered form. Furthermore, the 
observations of distant supernovae, at cosmological distance. 
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provide a direct evidence for an accelerating universe, which is 
naturally explained by the gravitational domination of a com- 
ponent with a relatively large negative pressure, Pq - WqPq 
with Wq < -1/3. The cosmological constant A (for which 
wa = -1) is historically the first possibility which has been 
introduced and which satisfies this requirement. However, the 
presence of a non-zero cosmological constant is a huge prob- 
lem in physics: (i) quantum field theory predicts that A should 
be the sum of a number of enormous contributions, so in order 
to avoid a cosmological catastrophe, it is usually assumed that 
a yet unknown mechanism produces a cancelation between all 
these contributions; (ii) it is difficult to think of a mechanism 
which puts A to exactly, but it is even more difficult to find 
a mechanism which gives pA ~ pc ~ as the super- 

novae observations suggest, where pp, and p^ are the Planck en- 
ergy density and the critical density today, respectively. For this 
reason the concept of quintessence, a scalar field with negative 
pressure, has recently been proposed as a possible alternative 
to a cosmological constant. 
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In this paper we shortly describe the quintessence paradigm 
and its effect on some observable quantities. We then sum- 
marise the different sets of data and methods used to constrain 
cosmological parameters. We then conclude by showing the re- 
sults on quintessence and cosmological parameters. 

2. Quintessence 

The idea of quintessence was proposed in order to allow the 
presence of an non zero dark energy, as suggested by obser- 
vations, without being confronted to the dramatic fine-tuning 
problem of the cosmological constant. Indeed, one still as- 
sumes that an unknown mechanism puts the bare cosmological 
constant is zero, and the smallness of the dark energy has a dy- 
namical origin, coming from a scalar field (p which has not yet 
reached the minimum of its potential V((p). 

A large number of quintessence models were already 
considered in the literature. Historically, the first proposed 
quintessence model had an inverse power la w pote ntial, V{(p) - 
M^^"_[f' (Ratra & Peebles |1988| ; Wetterich |1988| ; Caldwell et 
al. 1998), where the exponent a is positive and M is an en- 



ergy scale fixed so that the scalar field has the correct energy 
density today. This unusual shape, which can have some moti- 
vations from particle physics (Binetruy 1999 2000), insures 
that regardless of the initial conditions, the field will reach 
a so-called "tracking regime" (Steinhardt et al. 1999), where 



its pressure and energy density tend to a constant ratio given 
by Pq/Pq = Wq - (oWb - 2)/(a + 2), where Wb represent 
the pressure to energy density ratio of the other background 
matter fluids (photons, neutrinos, baryons, and cold dark mat- 
ter). Moreover, since we have in this regime Wq < Ws, this 
ensures that the quintessence energy density decreases more 
slowly than that of the background fluids and that ultimately, 
the quintessence field will become dominant. When this oc- 
curs, i.e., when its density parameter reaches Qq > 0.5, the 
field slows down in its potential, and reaches asymptotically 
Wq - -I (Steinhardt et al. 1999). The rate at which one goes 
from the tracking regime to the cosmological regime is usually 
quite slow, so that unless Wq is already close to -1 in the track- 
ing regime, it will still be significantly different from -1 when 
Q.Q ~ 0.7. For example, if cc = 6, in which case Wq = -0.25 in 
the tracking regime during the matter-dominated era, one has 
Wq - -0.4 today if Oq = 0.7. This feature leaves the hope of 
distinguishing a quintessence field from a cosmological con- 
stant. 

The effect of a quintessence field on CMB anisotropics is 
twofold (Brax et al. 2000). First, when the quintessence field 
becomes dominant it modifies the expansion rate of the uni- 
verse. This translates into a modification of the angular distance 
vs. redshift relation, and hence a shift in the peak structure of 
the CMB anisotropics power spectrum, the Cf's, for ( > 100. 
Second, the gravitational potentials decay at late time as the 
universe is no longer matter-dominated. This produces a so- 
called integrated Sachs-Wolfe effect and modifies the Cr spec- 
trum at low multipoles < 20) as a consequence of the fact that 
photons exchange energy with time-varying gravitational po- 
tentials. Both of these effects are also present with a cosmolog- 
ical constant, but they differ quantitatively with a quintessence 



field: the shift in the peak position is smaller, whereas the inte- 
grated Sachs-Wolfe effect can be very different (Caldwell et al. 



1998) 



Most of the CMB experiments do not cover a large fraction 
of the sky. On the co ntrary, t he new Archeops data (Benoit et al. 
2002a| ; Benoit et al. 2002b| ) are extremely precise around { ~ 
200 improving by a factor of two the precision measurements 
on the location of the first Doppler peak (Benoit et al. 2002c), 
whereas at larger angular scales (low {) the COBE data are 



limited by a large cosmic variance (Tegmark 1996). Therefore 
we can hope being able to contrain the quintessence parameters 
by through their influence on the position of the first Doppler 
peak rather than through the integrated Sachs-Wolfe effect. It 
is well-known that the position of the first peak is primarily 
sensitive to the curvature but also to several other cosmological 
parameters. It is therefore important to investigate a wide space 
of parameters in order to constraint the possible existence of 
quintessence in a robust way. 

Finally, let us add that quintessence also modifies sig- 
nificantly the matter power spectrum: as for the cosmologi- 
cal constant, matter fluctuations stop growing at the onset of 
quintessence domination. This has an influence on the normal- 
ization of the matter power spectrum on small scales, cr^. This 
effect of Wq can be understood as follows: as long as Wq is not 
too close to 0, one can roughly consider that the universe has 
experienced two distinct epochs since recombination: a first 
(usual) one where it was matter-dominated, and a second one 
where it is quintessence-dominated. The transition occurs when 
i^mat = fig = 0.5. If wc supposc that D.Q = 0.7 today, then the 
transition epoch occured at a redshift of 



,mat I 



- 1. 



(1) 



For Wq = -1 (cosmological constant case), one has z™' = 0.33, 
whereas for Wq = -0.2, this gives z™' = 3.1. If we suppose that 
perturbations grow as the scale factor before the transition and 
stop growing immediately afterward, then, this corresponds to 
a reduction of 



R 



1 + z™"' 



(2) 



in the matter power spectrum as compared to the cosmo- 
logical constant case. In practice, several factors modify this 
amount, such as the fact the the matter perturbation still grow 
(although slowly) in the quintessence-dominated era, but this 
illustrate the dramatic change in the matter power spectrum 
that quintessence can trigger [see also Benabed & Bernardeau 
(2001)]. When Wq - 0, the situation is even more dramatic as 



matter never dominates. In order to understand the value of erg, 
one is has to compute the growth of dark matter perturbations 
when one has a mixture of matter and quintessence with iden- 
tical equation of state. Neglecting the quintessence fluctuations 
(which are never subject to any instability), the Jeans equation 
is rewritten as 



'6 + 'H6- -'H^Q.^uatS = 0, 



(3) 
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Fig. 1. The amplitude of matter fluctuations erg as a func- 
tion of a constant Wq in cosmological models with Q.q - 0.7, 
Qb/i^ = 0.02, Hq = 66.5, n = I, and 0,o, 
normalized to match the COBE data. 



1, with amplitude 



When 6 represents the matter perturbation density constrast, 
an overdot denotes a derivative with respect to conformal time 
(defined hy adrj = df, where a is the scale factor and t is the 
cosmic time), and "H = a/a. The only difference with the usual 
matter-dominated case is that we do not have Qmat = 1 in the 
above equation. In the matter-dominated era, a oc rf^, so that 
1-1 -21 T], and one has 

5+ -6- ^^^,5 = 0. (4) 

When Qmat = 1, one gets the usual result that the growing 
mode goes as 6^ oc a. If we now suppose that Q^at = 0.3, 
then the growing mode goes as 6^ oc a" '**^. If this regime lasts 
for several orders of magnitude in redshift, then this represents 
an enormous reduction in the matter power spectrum, of around 
(1 H-Az)'-'", where Az represent the duration of this regime. This 
effect will imply that erg is drastically reduced when one goes 
from Wq = - 1 to Wq = 0. In order to illustrate this point we have 
computed the amplitude of matter fluctuations erg for a family 
of models which have the same same cosmological parameters 
but the equation of state Wq (Fig. [I]). As one can see erg varies 
dramatically with Wq and is therefore a very useful quantity to 
take into account in order to constrain quintessence. 

3. Methods and Data 

In the following, we make use of the most recent data avail- 
able on the CMB as well as on other relevant cosmological 
quantities in order to examine constraints that can be set on 
the amount of quintessence present in the universe. We assume 



Gaussian adiabatic fluctuations and a flat geometry of space. 
We assume a vanishing amount of gravitational waves as such 
a contribution has little impact on the position of the peak, but 
rather modifies the relative amplitude power between low and 
large £. Identically, a reionisation effect and a possible hot dark 
matter component are neglected in the following. We make the 
assumption that Wq - Const throughout all the epochs of in- 
terest. This assumption is unjustified in realistic quintessence 
models as one expects Wq to have varied at the radiation to mat- 
ter transition and to be decreasing today. However, if we want 
to study the influence of the quintessence field on the position 
of the Doppler peaks of the CMB anisotropies, then this is suf- 
ficent, and the constant Wq has to be seen as some "average" of 
a dynamical Wq(z) arising in a quintessence scenario. 



3.1. Likelihood from CMB 

In order to use CMB data, we first reconstruct the likelihood 
function of the various experiments. 

We follow the technique developped in Bartlett et al. 



(200C), and used in Douspis et al. (2001) and Benoit et al. 
(2002c), by constructing a large Cf power spectra database. We 
investigate six cosmological parameters assuming flat cosmol- 
ogy (Qtot - !)• The density of the universe is parametrised 
by the baryon contribution, Qb/i", and the dark energy, Qq, for 
which we allow its pressure-energy density ratio, Wq to vary. 

The Hubble parameter, Hq, the spectral index, n, and the 
normalisation of the spectra, parametrized in this work by erg 
are the remaining free parameters. Table |l] describes the corre- 
sponding gridding used for the database. 



We proceed by estimating cosmological parameters from 
the likelihood functions reconstructed as described in Benoit et 



al. ( [2002c[ ). We compute the value of the likelihood considering 
the actual band powers dataset of the COBE, BOOMERanG, 
DASI , MAXIMA, VSA, CBI, Archeops experim ents (T egmark 
Netterfield et al. |2002| ; Halverson et al |2002| ; Lee et 
Scott et al. I2OO2I; Pearson et al.l2002l Benoit et al. 



1996 



al. 2001 



2002b) on each model of our grid. In our approach, the best 
model is estimated as being those for which the likelihood is 
maximal Xmax, while the 68% (resp. 95%) 1 -parameter interval 
corresponds to -2 ln(£/Xinax) e [0, 1] (resp. -2ln(X/-Cmax) 6 
[0,4]) and the 68% (resp. 95%) 2-parameters interval corre- 
sponds to -2 ln(£/£inax) e [0,2.3] (resp. -2 ln(X/Xniax) e 
[0,6.18]). The likelihood shown in the following are already 
margilalized (by maximisation) over the calibration uncertain- 
ties and the amplitude. The results are presented as 2-D contour 
plots, showing in shades of blue the regions where the likeli- 
hood function for a combination of any two parameters drops 
to the levels corresponding to 68%, 95%, and 99% confidence 
regions. They would correspond to 1, 2, 3 er, respectively if the 
likelihood function were Gaussian. Dashed red contours mark 
the limits to be projected if confidence intervals are sought for 
any one of the parameters. To calculate either 1- or 2-D confi- 
dence intervals, the likelihood function is maximized over the 
remaining parameters. 
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Wq 




100 Qbh' 


Ho 


n 


o-g 


Min. 


-1.0 


0.0 


0.915 


25 


0.750 


0.1 


Max. 


-0.1 


1.0 


3.47 


101 


1.25 


1.3 


Step 


0.1 


0.1 


0.366 


xl.15 


0.015 


0.022 



Table 1. Grid of explored cosmological parameters; for Hq we 
adopt a log-periodic binning, Ho(i -H 1) = Ho(i) x 1.15 



3.2. Other Data 

In order to combine with other data of cosmological relevance, 
the corresponding additional likelihood have to be evaluated. 

An interesting useful additional constraint to add is those 
obtain on the amplitude of matter fluctuations on small scales: 
present day clusters data allow to constrain o-,., related to erg 
with a high accuracy, of the order of 5% (Blanchard et al. 
2000). Similar constraint can be obtained from weak lens- 



ing m easurements (Bacon et al. 200C; van Waerbeke et al. 
2002). However, significant differences among similar analy- 
ses have appeared in recent works based on clusters as well 
as on weak lensing measurements (Jarvis et al. 2002 ; Brown 
et al. 2002 , Hamana et al. 2002). We have therefore chosen 
to use two recent constraints, the differences of which will 
allow us to investigate a realistic range of systematic uncer- 
tainties. We first consider a constraint leading to high values 
of cTg (high cTg hereafter), in agreement with Pierpaoli et al. 
(pOOTl): o-s^",^. = 0-5 ± 10% (68% C.L.). Then, new estima- 



tions from Seljak et al. ( |2001| ), Viana et al. ( [2002| ), Reiprich 
et al. ( 2002 ), leading to lower values, are considered by taking 



o-sQ^^ = 0.43 + 10% (68% C.L.), hereafter low erg. The latter 
error bars estimates encompass the three above low normaliza- 
tion measurements of erg . 

We also use the recent determination of the Hubble pa- 
rameter from the Hubble Space Telescope (HST) Key Project 
(Freedman et al. [200 //q = 72 + 8 (68% C.L.), assuming 
Gaussian uncertainty. 

For the distant supernovae constraints we compute the like- 
lihood as follows. The magnitude-redshift relation of the super- 
novae is given by the following relation 



M = M + 5 logio £)z.(z, wg, Qq), 



(5) 



where M is the observed magnitude and Dl and Ai are respec- 
tively the "Hubble-parameter-free" luminosity-distance and the 
"Hubble-parameter-free" absolute magnitude at the maximum 
of the supernova. Dl is a function of redshift of the supernova, 
z, and of the cosmological parameters Wq and Qq. In the case 
of a flat universe with Wq as the pressure to energy density ratio 
of the quintessence field, the Dl function is given by 



£)l(z,Wq,Qq) = 
(l+z) 



'I 



du 



V(l - Qq)(1 + + Qq(1 + m)3(i+'^'q) 



(6) 



With the sample of supernovae from Perlmutter et al. ( 1999 ), 
we construct the likelihood Xsn with the following procedure: 



XsN = C exp 



n- 



(M*^-M,(Mwq,Qq)) 



2\ 



2cr? 



(7) 
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Fig. 2. Present CMB dataset likelihood contours in the 
quintessence paradigm. The sharpness of contours at Qq - 0.9 
is due to grid effect. 



where C is an arbitrary constant. We take for measured mag- 
nitude M°^^ and uncertainty cr the B-band peak magnitude and 
the total uncertainty of the supernovae from Perlmutter et al. 



(199S). At each point of the two parameters grid, Wq and Qq 



is associated a likelihood value obtained by maximising it over 
the absolute magnitude Al. 



4. Results 

4.1. CMB alone 

Constraints given by the CMB on some of our investigated 
parameters are shown in Fig. ^. Considering only CMB con- 
straints leads to degeneracies between parameters. Fig. ||shows 
the case of two parameters, n and Qb/i^, which are not af- 
fected by the assumed equation of state of the dark energy. 
Their prefered values and error bars are n - 0.95 + 0.05 and 
Qb/i^ = 0.021 + 0.003 (68% C.L.). Using CMB alone leaves 
the 2-parameters space (Qq, Wq) almost unconstrained. Finally, 
2D diagrams Qq vs. //q or Wq vs. Hq show significant level 
of correlation, but with degeneracies. This is illustrated by the 
plot Qq vs. Ho shown in Fig. ^ In our analysis, we found that 
with the improvement of CMB data obtained by the addition 
of Archeops band powers reduces appreciably the contours of 
constraints on the quintessence parameters as well as on cos- 
mological parameters because of the position and the ampli- 
tude of the first acoustic peak are better determined, but still 
does not allow to break the degeneracies. 
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Fig. 3. Likelihood contours with combination of CMB and priors. From left to right, panels show (i) colored contours correspond- 
ing to the combination CMB + erg (high), wiht overplotted SN contours; (ii) colored contours corresponding to the combination 
CMB + cTg (low), with overplotted SN contours; (iii) colored contours corresponding to the combination CMB + erg (low) + 
HST, with overplotted SN contours; (iv) colored contours corresponding to the combination CMB + erg (low) + HST + SN. 



4.2. Adding Non CMB priors 

As cluster abundance observations lead to a strong constraint 
on the normaUzation of the matter power spectrum, and erg is 
rather sensitive to the changing the equation of state, it is natu- 
ral to expect that this constraint in combination with constraints 
from Ce will lead to tight constraint on quintessence scenarios. 
Fig. H (left panels) shows the combination of CMB data with 
the two (Te priors described previously. As one can see, con- 
straints on the amplitude of matter fluctuations on small scales 
has the potential to break the degeneracies between Qq and 
Wq. Only a band-shaped region of the plane is not excluded. 
Furthermore, combinations with different priors lead to differ- 
ent likelihood contours, due to the strong effect of the equa- 
tion of state on erg emphasized in Section 2. More specifically, 
we find that the CMB, combined with the high normalization 
leads to a prefered region which is marginally consistent with 
the constraints given by high redshift supernovae (overplotted 
black Unes in Fig. The "concordance model" (Qq = 0.7, 
Wq - -1) lies on the 99% C.L. contours, and defines the two 
dataset as inconsistent. The best model appears then to have 
Qq - 0.4, Wq = -0.5, which is itself outside the 99% confi- 
dence region of the supernovae constraints. The combination of 
CMB with low normalization leads on the other hand to likeli- 
hood contours in agreement with those of supernovae and HST 
key project determination of //q. 

Due to the form of the joint CMB + erg (high or low) 
contours, a combination with the HST constraints is expected 
to give stronger constraints on both Oq and Wq. The corre- 
sponding likelihood contours of Fig. ^ (middle right panel) 
show that quintessence is not prefered over classical cosmo- 
logical constant even if the degeneracy is not totally broken: 
Qq = 0.70!°;;^, Wq = -l+O'* (95% C.L.). 

In order to break the degeneracy, it is clearly necessary 
to consider the additional information on the angular distance 
coming from distant supernovae. Considering a flat cosmol- 
ogy, the information on the luminosity of the supernovae can 
be expressed in term of constraints on the dark energy den- 
sity and equation of state. Prefered value are consistent with 



a cosmological constant, and the likelihood contours are al- 
most perpendicular to those of CMB, as shown in Fig. |[ 
Combining all the priors finally allows to put strong constraints 
on both quintessence parameters (Fig. ||, rightmost panel): 
Qq = 0.70!;}J^, Wq = -1+0-25 (95% L ) and finally breaks 
the (Ho, Qq) degeneracy, see Fig. Q 

As a main result, it appears that the classical ACDM sce- 
nario is then conforted and given the priors we used there is no 
need for quintessence to reproduce the present data, although 
quintessence models with low Wq are still viable, and that good 
fits to the data can also be found for models with w„ < -I 



(Melchiorri et al. 2002) However, those correspond to rather 
unusual quintessence models. 

Going back to typical quintessence models, finding only 
low values of Wq raises some interesting points. If we consider 
a pure inverse power law potential, having Wq < -0.6 when 
Qq = 0.7 implies a low value of the exponent a (if the slope of 
the potential is too steep, the field does not stop rapidly when 
it starts dominating, and Wq is too large). This translates into 
an annoyingly low value of the energy scale M arising in the 
potential, since one has M ~ Mp,(pc/m4)1/<4+o') (gj-ax et al. 
[2000 ). For example, in order to have M > 10"^ TeV, one needs 



a > 6, whereas data favor Wq < -0.6, which corresponds to 
a < 3, or M < 20 GeV. There are of course several possibil- 
ities which allow to have a lower Wq in quintessence scenar- 
ios. An example to evade this problem is to use the SUGRA 



potential found by Brax & Martin (1999), which is given by 
V - M'^^"/(p" X exp(i(;*^/Mp|). The exponential coiTection 
comes from supergravity considerations and insures that the 
potential has a local minimum which is almost reached by the 
field today. For example, for Qq = 0.7, one has Wq - -0.8, 
almost independently of a. This illustrates the fact that if it is 
still difficult to distinguish some quintessence models from a 
cosmological constant with the present state of available data. 



During the preparation of this work, an other g roup submitted a re 



lated paper to the archive (Melchiorri et al. 2003). Although the anal 



ysis is not exactly the same, their conclusions are similar. 
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Fig. 4. Likelihood contours with CMB + all priors. 



although these already allow to exclude a large number of 
quintessence models, among which the simplest one. 



5. Conclusion 

We have studied the constraints that can be obtained on cosmo- 
logical parameters within the quintessence paradigm by using 
various combinations of observational data set. For simplicity, 
only models with constant Wq were examined as they are be- 
lieved to be sufficient to existing data at their present accuracy 
level. For similar reasons, we neglected possible reionization or 
non-zero gravitational wave contribution. Indeed, this approach 
has the advantage to help to understand whether quintessence 
models should be favored or not over a classical cosmological 
constant, while constraints on more specific scenarios have to 
be investigated directly one by one (Douspis et al. 2002 1. Our 
analysis method has been to investigate contours in 2-D pa- 
rameters space. Such an approach allows to examine possible 
degeneracies among parameters which are not easy to identify 
when constraints are formulated in term of a single parameter. 
For instance, we found that CMB data alone, despite the high 
precision data obtained by Archeops do not require the exis- 
tence of a non-zero contribution of quintessence, because of 
the degeneracy with the Hubble parameter: in practice CMB 
data leave a large fraction of the (Qq, Wq) plane unconstrained, 
while only a restricted region of the (Qq, Hq) plane is allowed. 
On the contrary, we found that almost no correlation exist with 
the baryonic content Qb nor the primordial index n. In order 
to restrict the parameter space of allowed models, we have ap- 
plied several different constraints. Interestingly, we found that 
the amplitude of the dark matter fluctuations, as measured by 
clusters abundance or large scale weak lensing data can po- 
tentially help to break existing degeneracies, although existing 
uncertainties, mainly systematics in nature, do not allow firm 
conclusion yet. Clearly this will be an important check of con- 
sistency in the future. We have then added constraints from su- 
pernovae data as well as HST estimation of the Hubble param- 
eter in order to break existing degeneracies. This allows us to 
infer very tight constraints on the possible range of equation of 
state of the dark energy. Probably the most remarkable result is 
that no preference for quintessence does emerge from existing 
CMB data although, accurate measurement of the amplitude 



of matter fluctuations on scale of 8 /i ' Mpc may change this 
picture. 
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